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Abstract. We propose and analyse a Raman spectroscopy technique for probing the 
properties of quantum degenerate bosons in the ground band of an optical lattice. Our 
formalism describes excitations to higher vibrational bands and is valid for deep lattices where 
a tight-binding approach can be applied to the describe the initial state of the system. In 
sufficiently deep lattices, localized states in higher vibrational bands play an important role in 
the system response, and shifts in resonant frequency of excitation are sensitive to the number 
of particles per site. We present numerical results of this formalism applied to the case of a 
uniform lattice deep in the Mott insulator regime. 



1. Introduction 

The observation of the Mott-insulator state of a degenerate Bose gas in an optical lattice HI 
has opened an exciting new avenue for investigating strongly correlated condensed matter 
systems. Ensuing experiments (e.g. (2] [3] HJ) have investigated a wide variety of phenomena 
making optical lattices one of the leading systems for studying quantum atom optics. 

In the initial experimental report, two key pieces of evidence were given for the Mott 
insulator state: the loss of phase coherence, and the appearance of a gap in the excitation 
spectrum. More recent experimental work has used Bragg spectroscopy ||5] to study the 
equilibrium properties of this system. While this type of probing reveals the appearance of 
a gap in the excitation spectrum, experimental applications of this technique operate well- 
beyond the linear response regime, making direct comparison with theory difficult. 

In this paper we propose and theoretically analyse a Raman spectroscopy scheme for 
probing Mott insulator states in the lattice system. The formalism we present is sufficiently 
general to include the excitation of atoms to higher vibrational bands where they should 
be easily discernible from the unscattered atoms in time-of-flight analysis. We apply this 
formalism to a uniform lattice and show how the correlated nature of the Mott insulator state 
can give rise to localized states in the excited bands when the lattice is sufficiently deep, and 
that the resonant frequency of these localized states is sensitive to the local filling factor in the 
optical lattice. 

| The observable for Bragg spectroscopy in optical lattices is the energy transferred to the system. This observable 
is rather difficult to measure accurately, and to obtain a signal experiments necessarily add a large amount of energy, 
well-beyond the linear response limit. 
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Figure 1. Schematic diagram of the spectroscopic process. Atoms in the ground band of 
system 1 are transfered into the initially unoccupied system 2 by a two-photon Raman process 
through a virtual level \v). A residual energy difference liwo and momentum kick tu\ (not 
shown) are transfered to the atoms by each scattering event. 



2. Formalism 

The system of interest, which we will refer to as system 1 , is a degenerate collection of bosonic 
atoms populating the lowest vibrational band of an optical lattice well-characterized by the 
Bose-Hubbard Hamiltonian 

Hi = e ™j ■ - J H 4 & j + Ull ® w - !)> (!) 

where ctj is a bosonic operator that annihilates an atom from the Wannier state Wj (x) centered 
on lattice site j, with fij = fiteij the respective number operator. The quantity J, known as the 
tunneling matrix element, characterizes the tunneling between lattice sites and is determined 
from band structure calculations (7][8). Interactions between particles are described by the 
matrix element a w = J d 3 x |w,,(x)| 4 and the coefficient Uu — Airanh 2 /m, where an is 
the s-wave scattering length for collisions between atoms in internal state 1. We restrict our 
attention here to the translationally invariant system (i.e. neglect the influence of external 
trapping potentials in addition to the optical lattice), and the constant eo characterizes the 
mid-point energy of the ground band. 

In this work, we consider a scheme for probing the properties of system 1, by an internal 
state changing Raman transition, of the form used in Ref. [9| to output couple an atom laser. 
A theory for this type of process in an optical lattices has also been proposed by Konabe 
et al. [25 (f| This is in contrast to a recent theory iflOl [TT1 [121 [131 [T4l [TBI and experiments 
|EJ [6] that have considered internal state preserving transitions to perform spectroscopy - 
known as Bragg scpectroscopy [16). A schematic diagram of this coupling is shown in Fig. 

§ We note that the primary difference of our treatment to that of Konabe et al. , is that we include interactions between 
scattered and unscattered atoms, which are essential for the main results we present here. 
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[T] Following the formalism in Ref. IfTTl this type of coupling is described by an interaction 
term of the form 

V = YpQ( t ) I d 3 xV»t(x)a J w J (x)e i ( q - x -" ot ) + H.c, (2) 
2 . J 

where is the field operator for bosonic atoms in the second internal state, and the quantities 
h<\ and hujQ, specify the momentum and energy transfer of the Raman process respectively. 
For simplicity we define Hluq to be the excess energy transferred over the internal state energy 
difference. We take the amplitude of the Raman coupling to be of strength V p beginning 
at t = 0, and of duration T p . We will assume that the internal states 1 and 2 are different 
hyperfine ground states for which we can neglect any collisional spin evolution (see |[T8l[T9l ). 
To arrange such a Raman probe experimentally, a pair of light fields (in addition to those used 
to create the lattice) with appropriately chosen polarizations to couple the states of interest 
will need to be applied to the atoms in the lattice. While such a probe could be focused down 
to address a few sites in a lattice, our theoretical development here is for the case where the 
Raman fields uniformly illuminate the system and have no spacial selectivity. 

Our interest lies in the linear response regime, where only a small portion of the atoms 
are scattered into internal state 2. The evolution of these atoms in internal state 2, which we 
will refer to as system 2, is described by the Hamiltonian 

H 2 = J d 3 x^|(x)(tf sp + U 12 %-K(x)| 2 )V> 2 (x), (3) 

3 

where 

H sp =p 2 /2m + V cxt (x), (4) 

is the single particle Hamiltonian with V ex t (x) the external potential (taken to be the same 
lattice potential experienced by atoms in system 1), and U12 = Aira^H 2 /m characterizes 
the interactions between particles in systems 1 and 2. Within the validity regime of linear 
response it is permissible to neglect interactions between atoms in system 2, as their density 
will remain low. 

System 2 is initially in the vacuum state, and the spectroscopic signal to be measured is 
the number of atoms in system 2 after the Raman pulse. This is a convenient experimental 
observable as detection techniques can readily distinguish between atoms in different internal 
states [ 1 9 1 . Additionally when the coupling produces atoms in an excited vibrational band, the 
excited atoms can be differentiated by their behavior upon expansion from the lattice lEOlETI . 
These properties of the observable should enable the system response to be measured for small 
amounts of Raman excitation, allowing the system to be probed in the linear response regime. 

The expression for the number of excited atoms, derived using linear response theory, is 

^2 > = " / duR(q,uj) — -£ , (5) 

where we have introduced the correlation function 



x J>*(xH(x') (aJ(t)^ 2 (x,t)^( x ',o)a i (o) 



(6) 



Eq 
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In this expression the time dependence of the operators is in an interaction picture with respect 
to unperturbed Hamiltonians Hi + H 2 , and the expectation is taken on the initial equilibrium 
ensemble. 

We will now show how to evaluate the correlation function i?(q, to), and how it relates to 
the properties of the atoms in system 1. We take the initial density matrix to be p = p\ ® p™ c , 
where pi is the degenerate lattice state we are interested in probing and p™ c is the initial 
vacuum state for system 2. We can then show that without further approximation, i?(q, uj) 
can be written in the form 



R(q,u) = l-jdtjd 3 X J 



d 3 x / el q-(x'-x) 



J2(q aJ(*)%(o)Pi q)(o 2 ^ Q (x,i)v;|(x',o) o 2 )<(xh- 

Qij 

where the variable Q represents a trace carried out over the number state basis for the operators 
{a,}, i.e. Q <-> | . . . , n®_ 1 ,nf , . . .), and IO2) represents the vacuum state of system 2. The 
operator $9 (x, t) is -02 (x) evaluated in an interaction picture with respect to the Hamiltonian 

Hf = f d 3 x^(x)(iJ sp + U u ^2(nf - (5y)| Wi (x)| 2 )v3 2 (x). (8) 

1 

We note that in deriving Eq. (0 we have made use of the result 

e i ^*/^ 2 ( x ) e -^ 2 t/ft ai(0) | Q) = a 3 ( )|g)V,Q(x,i), (9) 

where the 5ij term in Eq. ([8} arises from commuting cij with exp(iH 2 t/h), and we have made 
the replacement hi — > nf as \Q) are number states. 

We refer to as the Q-defect Hamiltonian for system 2, as it arises from the removal 
of a system 1 atom from site j of the number state Q. As the system 1 atoms form an effective 
potential for those in system 2, the removal of an atom at site j due to the Raman excitation 
creates a potential hole (i.e. the nf — 5ij term in Eq. ([H}). This defect plays a key role in the 
response spectrum as we will demonstrate later. Since Hj is a quadratic Hamiltonian it can 
be diagonalized by numerical methods to obtain its eigenvectors {^^(x)} and eigenvalues 
{hujf m }, i.e. fiuf m (j>j m (x) = H jQ (t)f m (x), where m is the quantum number specifying the 
state and j labels the defect location. We note that in the limit U12 — » 0, the cj)^ m reduce 
to the Bloch states of H sp and m becomes the quasimomentum and band index. We obtain 

ipf(x, t) — J2 m ^fm( x )^L e ^ mt > wnere b^ m is a bosonic annihilation operator and arrive 
at the expression 

i?(q,c) = -L J dte^^T.of^A^A^e-^, (10) 

Q ijrn 

where we have defined 

A% = J d 3 x^tWe 4q - x ^(x), (11) 

cg(t) = (Q\al(t)aj(0) Pl \Q). (12) 

Eqs. d2J and ([Tol l represent the key results of this work, and we now briefly comment 
on the physical process they describe. Fundamentally, i?(q, u>) characterizes the excitation 
spectrum of atoms from system 1 into system 2. In the context of ultra-cold gases, a result 
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similar to our starting point [Eq. Q] has been given by Luxat et al. Il22l for the case 
of a harmonically trapped Bose gas (also see 11231 12410 . However their treatment neglects 
interactions between atoms in different hyperfine states and assumes the single particle 
correlation functions for the atoms in systems 1 and 2 are independent. The extension of 
this theory to the optical lattice has been provided by Konabe et al. Il25l . For the current 
experiments with Rubidium atoms in a deep optical lattice these approximations are not 
tenable. In particular, as we noted previously, the lattice site from which the atom is excited 
acts as the localizing defect and accounting for correlations between the systems (as we have 
done in Eq. (O) is essential. 

We note that our formalism [i.e. Eqs. © and (TT~0b l is quite general. As long as the 
dominant number states of the many-body state are known, e.g. through exact diagonalization 
or Matrix Product Decomposition techniques (e.g. esee 11261 ). then the Raman response can 
be determined. In the following sections we consider the application of our formalism for 
two special cases. The results we present in the next section are calculated for ID systems 
for numerical convenience, though our interest is in the regime where the scattering between 
particles is three-dimensional and well-described by a contact interaction. 

3. Single Site Limit 

The limiting case of a single tightly confining harmonic well of frequency Who is a useful 
approximation to a deep lattice in the regime where tunneling between sites can be neglected. 
This limit shows the main physical features of Raman spectroscopy and provides a useful 
approximation to the full solution. Assuming that interaction shifts are small compared to 
the oscillator energy fio>ho, we may approximate the modes of the system as being harmonic 
oscillator eigenstates {(/? m (x)}, with respective energies {e m — frui^^m + 1/2)}. For the 



single site case our formalism maps according to 

Wj(x) ->-^ (x), (13) 

-> *rn{*), (14) 

A% ^A m E|d 3 x^(x)^> (x), (15) 

ttuf m -> Tiw™ = e m + U 12 a Qm (n - 1) (16) 



where ao m = f d 3 x\ipo(x.)\ 2 \ip m (x.)\ 2 . Since the Hamiltonian for system 1 reduces to 
Hi — > eoh + Uiia w n(n — l)/2 in this limit, we have replaced Q by n (i.e. the single 
site many-body number states are just single mode number states n}) and evaluated the 
correlation function as 

cg(*) -> <&•(<) = «exp(i[e + U n a w (n - l)]t/h). (17) 
Using Eq. ([TOj with p x = ^ nn , \n)p nn ,(n'\, we obtain 

R(q.,u) = ^2p nn n\A m \ 2 S(uj ICS (n) - ui), (18) 

71 

hhj rcs {n) = e m - e + [Ui 2 a 0m - Una w ]{n - 1). (19) 

We have also assumed that the probe only couples to a particular excited state m, with the 
other states sufficiently far detuned that their contribution is negligible. 

Eq. ( fT8l shows that the response of the system is proportional to p nn , and most notably, 
if [Unworn — Ui\a w ] 7^ 0, then the response frequency is linearly dependent on the value of n 
(i.e. the number of atoms at the site). In this regime the Raman spectrum reveals the number 
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distribution at the site. For the case of 87 Rb (the atom of primary interest in bosonic optical 
lattice experiments) the interactions between the relevant hyperfine states are approximately 
degenerate (i.e. U\\ ~ E/12) > so that the difference a 0m — a w will be the primary factor in 
determining the magnitude of the number dependent shift. As an immediate consequence we 
note that for the case U\\ = Un and a Raman pulse coupling to ground vibrational state of 
system 2, then the spectrum will be independent of n (i.e. for m = we have ao™ = a w ). 
Therefore to obtain a number dependent spectral response with 87 Rb will require scattering 
into excited vibrational states of system 2. 

4. Uniform Mott Insulator 

We now consider the more general case of a translationally invariant lattice with N s sites 
and periodic boundary conditions. We assume that the number of atoms in the system is 
commensurate with the number of lattice sites and system is deep in the Mott insulating 
regime, where Uua w 3> J. In this regime the many -body ground state is well approximated 
as \Qn) = \ • ■ ■ ,n,n, . . .) (i.e a definite number of atoms at each site). In calculating 
the response of the system to the Raman probe according to Eq. ©, the summation 
over Q reduces to this single state. The next order correction to the translationally Mott 
state is particle hole states fl3l . which contribute to the ground state with an amplitude 
J/Una w <C 1, and can be neglected. 

To obtain an approximation for the temporal correlation function of system 1 in state 
\Q n ) we ignore the tunneling term in H\. This approximation amounts to neglecting particle 
tunneling between sites in system 1 over the time scale of the Raman probe, and should be a 
good approximation in deep lattices. In this limit the correlation function is 

c9f(t) = n exp(i[e + U u a w (n - l)]t/h)S. lv (20) 

and making use of the translational invariance in evaluating Eq. dTOT t we obtain 

R( q ,w) = £nJV.|^| a %-a£r), (21) 

tn 

"?es m = ~ [ £ o + U n a w (n - l)]/h. (22) 

Note that due to translational invariance the precise value of I used in Eq. (|2"TT i is unimportant. 

We now calculate the response spectrum {^2^ using Eqs. © and (|2TT > for a translationally 

invariant ID lattice. We take the lattice potential to be V eK t{x) = Vd cos 2 (x/\), arising from 
counter-propagating lasers fields of wavelength A, and we specify the lattice depth Vd in units 
of En — h 2 /2mX 2 . The calculation is performed for typical 87 Rb parameters (e.g. see Ref. 
[T)) to demonstrate the practicality of our probing scheme. We will not consider the influence 
of q on the spectrum in this papei|jjl and for simplicity we have taken q as being identical to a 
reciprocal lattice for our calculations. 

The results we have obtained are shown in Fig. [2] For the case of Vd = 10-Er (Fig. EJa)) 
the Raman response is shown over a frequency range that includes resonant coupling to the 
first two excited bands. The superimposed graphs show the spectra for Mott insulating states 
of various filling factors and clearly exhibit frequency shifts proportional to n. Additionally, 
we notice that in the first excited band (12kHz - 18kHz), the dominant response peak occurs at 
the low end of the spectral feature (e.g. the peak at ~ 15.5kHz for n — 1), adjacent to a broad 
base. This is feature is also seen in the 2nd excited band (24kHz - 36kHz), however the peaked 

|| Generally large q favors coupling to higher bands, though certain choices of q can suppress coupling efficiency, 
e.g for q = the response from the first band is zero by symmetry. 
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Figure 2. The percentage of atoms excited after Raman excitation in a ID lattice of depth (a) 
Vb = 10-Er (b) Vj3 = 30Er. In each plot the spectrum of Mott insulator states with filling 
factors of n = 1 (dark line), n = 2 (medium line) and ra = 3 (light line) are shown. The 
respective single site results for m = 1 are given as dashed lines. Parameters are q = tt/X, 



T p = 1.5ms, V p 
harmonic with fx 



0.05E R , A = 850nm, N s 
! 37.6kHz, and an = ai2 = 



= 51, transverse confinement taken to be 
5.29nm. 



feature is much less dominant relative to the broad base. The peak originates from an excited 
band state that is partially localized above the defect site. This localization leads to a strong 
coupling matrix element and as this state resides at the defect, its energy is lower than 
the other states in the excited band. We refer to this state as the defect state, which is clearly 
identified as the most strongly excited feature at the bottom each spectral band. Since the 
localization is not perfect, the other states in the excited band have appreciable amplitude at 
the defect site. These states give rise to the broad though more weakly excited, band of states 
above the resonant peak. In the 2nd excited band, the same features are seen, however as the 
effective tunneling in this band is much higher, and the defect state is much less localized. 

In Fig. |2jb) the response spectrum is shown for a lattice of depth Vd = 30-Er. At this 
depth the response from the first excited band states have shifted up to 30kHz, and the second 
excited band has moved out of the frequency range considered. In contrast to the spectrum in 
Fig. El a), we only notice the resonant peak due to the defect state, without a discernible broad 
base of band states. This arises because at this depth the tunneling rate in the excited band 
is sufficiently small that the defect state a becomes completely localized. The other states in 
the excited band are necessarily orthogonal to the defect state and so have vanishing coupling 
matrix elements Az°,, . 

For comparison in Figs. |2ja) and (b) we also show the single site predictions for 
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the spectrum calculated using expression ( fT8l with c^ho chosen to match the effective trap 
frequency at the lattice site minima. The frequency location of the response spectra compares 
badly with the full lattice solution, arising primarily from the inadequacy of the harmonic 
approximation for accurately predicting band structure. The location of the spectra can be 
easily corrected for by calculating the term e m — eq in the expression for u) les (n) using the 
non-interacting band structure result. However, the single site approximation captures many 
of the salient features of the full lattice solution, such as the magnitude of the n-dependent 
shift in the the response spectrum. In the strongly-localized defect limit (Fig. EJb)), the single 
site approximation quantitatively predicts the response amplitude as the role of the band states 
can be neglected. 

To quantify the emergence of the defect state we examine the spectrum of w^" 1 as 
a function of Vd in Fig. Oa). For Vd ^ 13.5-Er a single defect state is observed to 
drop below the first excited band. The condition for the emergence of a strongly localized 
defect state is that the energy reduction from localization above the defect in system 1 is large 
compared to the effective inter-site tunneling in the excited band. To verify this criterion we 
compare the excited bandwidth (characterizing the excited band tunneling rate) and energy 
reduction at the defect, approximated by Uua w . These two energy scales are seen to cross at 
Vd ~ 13.5-Er. Finally, we note the validity condition for linear response treatment of Raman 
spectroscopy. This requires the number of atoms excited to remain small compared to the 
total number of atoms in the system. In terms of the Raman parameters this condition is given 
as nVpTp /4h 2 <C 1, where n is the average number of atoms per site. 
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We briefly comment on the relationship of the ID calculations presented here to an 
equivalent system in a 3D lattice. The primary difference in applying our formalism to a 
fully three-dimensional system is that Eq. ([8]) will need to be solved in 3D. For the case of 
coupling to excited bands, the tunneling between lattice sites can occur in all directions, but is 
dominated along the direction in which the vibrational excitation has occurred. This direction 
will be parallel to the direction of the momentum transfer in the Raman coupling, which we 
take to be parallel to a lattice vector. The tunneling in the orthogonal directions will be given 
by the ground band tunneling rate which is typically much smaller. This suggests that the 
additional shifts in a fully 3D lattice will be of order the ground state tunneling rate and will 
thus contribute small corrections to the results presented here. A full study in 3D will be the 
subject of a future investigation. 

5. Conclusions and outlook 

In this work we have proposed and analysed a Raman spectroscopy technique for probing the 
properties of quantum degenerate bosons in the ground band of an optical lattice. We have 
observed that for sufficiently deep lattices, localized states in higher vibrational bands play 
an important role in the system response, and shifts in resonant frequency of excitation are 
sensitive to the number of particles per site. While our main study has considered the case 
of a perfect Mott insulator in a translationally invariant lattice, our results suggest that in the 
limit of strongly localized defect states the response of the system is well-described by the 
single site result, and thus only depends on the local number distribution at each lattice site. 
The Raman spectrum may therefore be a useful method for measuring the relative portion of 
system 1 at sites with filling factor n atoms. We speculate that in the strongly localized limit, 
the homogeneity of the lattice is rather unimportant, with the existence of the defect states 
arises from the large difference in the effective potential energy between the site where the 
particle is excited and the neighbouring site. This suggests that many of the predictions we 
have made here, and in particular the results of the 1-site model, should qualitatively apply to 
inhomogeneous lattices, such as the combined harmonic and optical lattice potentials made 
in experiments. In this case, as well as in the superfluid limit (where significant number 
fluctuations exist), significant corrections may arise from resonances between neighbouring 
sites that would allow the defect to be localized over several sites. However, it seems 
reasonable to expect that these resonances would contribute to the broad background in the 
Raman spectrum, and the sharp features from Mott-insulating regions (if present) would be 
clearly visible. Characterizing the role of superfluid fluctuations and the external confining 
potential will be the subject of future work. 
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